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We consider doping of non-magnetic impurities in tlie spin-1/2, 1/5-depleted square lattice. This 
structure, whose undoped phase diagram offers both magnetically ordered and spin-liquid ground 
states, is realized physically in CaV409. Doping into the ordered phase results in a progressive 
loss of order, which becomes complete at the percolation threshold. By contrast, non-magnetic 
impurities introduced in the spin liquids create a phase of weak but long-ranged antiferromagnetic 
order coexisting with the gapped state. The latter may be viewed as a true order-by-disorder phe- 
nomenon. We study the phase diagram of the doped system by computing the static susceptibility 
and staggered magnetization using a stochastic series-expansion quantum Monte Carlo technique. 



PACS numbers; 75.10.Jm, 75.30.Hx, 75.50.Hm 

Doping of low-dimensional quantum spin liquids leads 
to a variety of phenomena, including unconventional su- 
perconductivity when localized spins are replaced by mo- 
bile holes Q|, and effective random spin systems, whose 
low-temperature properties are dominated by large ef- 
fective spin degrees of freedom, in the presence of 
static impurities [^j-Q. Such low-temperature phases 
have been observed recently in quasi-one-dimensional 
(Id) spin systems. In the two- leg ladder compound 
Cai4_a;Sra;Cu2404i, a superconducting transition has 
been measured at high pressure pj. Zn-doping of the 
quasi-ld ladder material SrCu203 [p|, and of the dimer- 
ized chain CuGeOa [Q, yields low-temperature phases 
with weak antiferromagnetic (AF) order, arising from ef- 
fective couplings of random strength and sign between 
unpaired spins. 

The question naturally arises of whether similar phe- 
nomena can be expected in quasi-2d systems. There 
are several prominent examples of 2d spin liquids, in- 
cluding (i) structurally frustrated materials, such as 
SrCu2(B03)2 1^, containing AF coupled spin-1/2 trian- 
gles, (ii) interaction-frustrated systems, such as the J1-J2 
Heisenberg model on a square lattice in the gapped 
regime around J2/J1 = 0.5, and (iii) structurally dis- 
cretized materials, such as the 1/5-depleted AF square 
lattice formed by the spin-1/2 vanadium ions of CaV409 
(Fig. 1) ||l^. These systems have resonant- valence-bond 
(RVB) ground states with a finite spin gap. Their AF 
correlation length is finite, decreasing with increasing 
gap size, and the low-temperature uniform susceptibility 
is exponentially activated. We find that non-magnetic 
impurities in such systems act to induce a long-ranged 
AF order, which exists on the subsystem of uncombined 
spins, and coexists with the dominant, gapped phase of 
the spin liquid. This order arising from disorder is the 
central result of our analysis. As in the quasi-ld case 
[0| , its origin lies in the unfrustrated nature of the ex- 



ponentially weak interactions between unpaired spins in- 
troduced when impurities break the RVB singlet units 
in a bipartite lattice: these are ferromagnetic when the 
free spins lie on the same sublattice, and AF for opposite 
sublattices 




FIG. 1. The 1/5-depleted square lattice. Intra-plaquette 
couplings J are indicated by solid lines, and inter-plaquette 
couplings J' by dashed lines. 

In this work we consider the effects of non-magnetic 
impurities in the 1/5-depleted square lattice geometry 
(Fig. 1), motivated by the structure of CaV409, and by 
the possibility of doping this particular material with Ti. 
The set of superexchange parameters required to model 
CaV409 has been found [^jl3| to be dominated by the 
frustrating, next-neighbor AF interaction, which to a rea- 
sonable approximation allows the real system to be con- 
sidered as an ensemble of weakly coupled metaplaquettes 
(a square of side \/2 larger than the nearest-neighbor 
spacing) of 4 V^^ spins. We present first a theoretical 
analysis for the unfrustrated system, which contains all 
the qualitative features we wish to address, and return 
thereafter to the experimental situation. 

The effective magnetic Hamiltonian, 

H = J^S, -S.+J'^S, -S,, (1) 

□ □-□ 

is a spin-1/2 AF Heisenberg model with J is the intra- 



1 



plaquette and J' the inter-plaquette exchange coupHng 
constant, in which the sums marked by □ and □ — □ run 
over nearest neighbors within and between the plaque- 
ttes. We win explore the full parameter range J'/ J > 0, 
which contains a spin-liquid, plaquette RVB (PRVB) 
phase for J'/J < 0.93 , an AF-ordered phase in the 
approximate interval J'/J G [0.93,1.68], and a further 
spin liquid, the dimer RVB (DRVB) phase, at larger 
J'/J > 1.68 Q. Expansions about the plaquette and 
dimer limits will be compared with simulations using the 
stochastic series-expansion quantum Monte Carlo (SSE 
QMC) algorithm |jl^, on 1/5-depleted, square lattices of 
up to 30x30 sites, and at temperatures down to T = 
0.02J. This algorithm involves expansion of the parti- 
tion function in inverse temperature, employs local and 
global system updates, and is significantly more efficient 
than conventional QMC schemes. In the presence of 
non-magnetic impurities, ensemble averages over typi- 
cally 500 random realizations were taken to ensure nu- 
merical errors smaller than the symbol sizes in Figs. 2 
and 3. 




FIG. 2. (a) Ensemble of decoupled plaquettes: tem- 
perature dependence of Tx{T) at impurity concentrations 
p = 0.0, 0.05, 0.4. Solid lines are exact results, and symbols 
are QMC simulations. Tx{T) is normalized by the number 
of spins in the system, (b) As (a), for decoupled dimers. In- 
sets: dependence of zero-temperature Curie constant C on p. 
Solid lines denote decoupled plaquette and dimer limits, and 
dashed lines the results from perturbation expansion with (a) 
J'/J = 0.1 and (b) J/ J' = 0.1. 

Considering first an ensemble of completely decoupled 
4-spin plaquettes, a small number of randomly placed, 
non-magnetic dopants creates free spin-1/2 moments, 
corresponding to the doubly degenerate — ±1/2 
ground state of 3-spin Heisenberg clusters Each 
of these "pruned" plaquettes gives a paramagnetic con- 
tribution to the susceptibility x(T) at low temperature. 
Upon further doping, clusters containing n < 3 spins 
are also created, with distributions P{n) |]l^, but only 
unpaired spins contribute to the low-T Curie tail. Fig. 
2(a) shows Tx{T) for various impurity concentrations. 
Without impurities, this quantity extrapolates to zero as 



the temperature is lowered, indicative of the spin gap in 
the clean system. With impurities, the zero-temperature 
Curie constant is given by 



C 



^hn[rx(T)] 



P(l) + 1/3P(2) 
1-p 



my 



(2) 



and depends on concentration as shown in the inset of 
Fig. 2(a). For all p, three distinct temperature regimes 
are observed, (i) For T ^ Amax (not shown) the spins 
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are essentially independent, each contributing to C . 
Amax — 1-5 J is the largest energy level separation in any 
cluster, (ii) For A„iin < T < A^axj some spins freeze 
into singlets, gradually reducing the overall magnetic re- 
sponse by a factor of approximately p as T decreases. 
Amin = 0.25 J is the smallest energy spacing in any clus- 
ter, (iii) For T < Amin, plateaus appear in Tx(T) when 
the contributing clusters reach their ground states. 
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FIG. 3. Temperature dependence of Tx{T) at impu- 
rity concentrations p — 0.0, 0.05, . . . , 0.4. (a) PRVB regime, 
J'/J = 0.3; (b) ordered AF regime, J'/J = 1.2; (c) DRVB 
regime, J'/J — 2.0. Solid lines are results from expansions 
about the decoupled limits, and symbols are from SSE QMC 
simulations. Dashed lines in (b) are guides to the eye. 

With finite inter-plaquette coupling, spin-1/2 moments 
on adjacent plaquettes may combine to form new RVB 
clusters with smaller spin gaps, involving larger numbers 
of spins. In the susceptibility, this effective recombina- 
tion reduces C from its value for decoupled, pruned pla- 
quettes, further suppressing Tx{T) at very low T. This 
cluster formation may be quantified from the distribu- 
tion functions for neighboring plaquettes, evaluated per- 
turbatively in J'/J. These functions have a non-trivial 
dependence on the relative positions of the dopant vacan- 
cies in adjacent clusters. Close to the PRVB limit, the 
inset of Fig. 2(a) illustrates this reduction of C, which 
is strongest at smaller p where more spins are available 
to form larger RVB clusters. Fig. 2(b) shows analogous 
results for the DRVB limit. 

In Fig. 3, the T-dependence of the uniform suscep- 
tibility is shown at various doping levels in the three 
regimes of J'/J. In Figs. 3(a) and (c), results from SSE 
QMC simulations (symbols) are plotted along with those 
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from perturbative expansions about the limits of com- 
pletely decoupled plaquettes and dimers (solid lines). At 
p = 0, Tx{T) clearly displays activated behavior in the 
two gapped regions, but a linear form in the long-range- 
ordered phase. 

PRVB: At small p, a plateau develops in Tx{T) at low 
temperatures, T ^ Amin, indicating a freezing of the pla- 
quettes into their ground state configurations. For a fi- 
nite inter-plaquette coupling, the additional reduction of 
C observed at ultra-low temperatures indicates that the 
plaquette spins become correlated on an ultra-low energy 
scale Jeff. At sufRciently small p, this energy can be esti- 
mated [H by assuming RKKY-like interactions between 
the impurity spins, yielding 

where S.af is the short-range AF correlation length in the 
spin-liquid regimes. Approaching T = 0, the Curie con- 
stant extrapolates to C = p/12, a reduction of roughly 
1/3 compared to the plateau value At larger impu- 
rity concentrations, the onset of this ultra-low-T regime 
moves up, consistent with T ~ Jeff (Eq. 3). At still 
higher doping, the picture of isolated impurity spins, and 
RKKY interactions mediated via the spin-liquid RVB 
background, breaks down. The regimes of moment for- 
mation and correlation overlap, and the plateau feature 
in Tx{T) becomes less pronounced. On doping towards 
the percolation threshold pc, the magnetic response is 
dominated by the largest percolation clusters, typically 
containing several pruned plaquettes. From simulations 
in the classical limit [|l8|, we find pc = 0.26(9) for the 
1/5-depleted lattice. Studies of percolation thresholds 
in regular square lattices jl^] suggest that the value of 
Pc for the quantum spin system differs by less than 1% 
from the classical value (pc = 0.407). In this regime the 
plaquette plateaus disappear, a trend accentuated as the 
long-range-ordered regime is approached {J' /J 0.93), 
where £^af increases and (Eq. 3) the onset of the ultra- 
low-T regime moves to lower energies. 

DRVB: As in the PRVB regime, non-magnetic impu- 
rities introduce local moments in the gapped, spin-liquid 
state. A dimer "plateau" is reached at T ~ Jog, and 
further moment recombination occurs at ultra-low tem- 
peratures, leading to an additional reduction of Tx{T) 
as T ^ 0. The plateaus disappear on approaching the 
transition to the AF phase, J'/ J — > 1.68+. 

AF: In contrast to gapped spin liquids, isolated, non- 
magnetic impurities do not introduce quasi-free local 
moments, because all spins remain aligned with the 
staggered AF background |^9|. However, the low- 
temperature susceptibility does have a divergent contri- 
bution, albeit with C rather smaller than in the spin liq- 
uids at given p, arising from free moments isolated from 
the AF system by dopants. 
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FIG. 4. Staggered magnetizafion as a function of impurity 
doping for coupling ratios J'/ J in the PRVB (dashed), or- 
dered (solid) and DRVB (dotted) regimes. The inset shows 
the evolution of with J'/J at a fixed doping p = 8%; 
arrows indicate the phase boundaries of the pure system. 

Returning to the issue of order by disorder, our simu- 
lations close to the AF regime of the pure system show 
that introduced moments do not simply recombine into 
larger RVB clusters, but form a long-range-ordered AF 
network at ultra-low temperatures. This doping-induced 
long-range order is best characterized by the staggered 
magnetization m^, which is shown in Fig. 4 as a func- 
tion of non-magnetic impurity doping within each of the 
distinct phases of the system. As also in Fig. 5, these 
data were obtained by finite-size extrapolation at ultra- 
low temperatures, and thus represent the thermodynamic 
limit. In the long-range-ordered AF regime, we see that 
random static vacancies simply reduce the average stag- 
gered magnetization, driving it to zero beyond the per- 
colation threshold pc — 0.26(9). By contrast, in the spin- 
liquid regimes the doping-induced moments have effective 
interactions mediated by the RVB background, which 
are unfrustrated A long-ranged order may then be 
expected to emerge from the large-spin clusters at zero 
temperature. Indeed, the resulting AF network has an 
extensive staggered magnetization, which peaks around 
p — 0.08. This ordered AF phase coexists [0 with the 
gapped, spin-liquid state of the majority, undoped pla- 
quettes or dimers (Fig. 3). The mechanism leading to 
this moment formation may be viewed as an order-by- 
disorder phenomenon. The inset of Fig. 4 shows the 
staggered magnetization as a function of the coupling 
ratio at the maximal doping. The disorder-induced, or- 
dered moments in the spin-liquid regions are significant 
over a wide range of J'/J, and become large near the 
phase boundaries of the pure system, where they cross 
continuously to the full moment of the intrinsically or- 
dered regime. 

Our results are summarized by the phase diagram in 
Fig. 5, where the shading illustrates the strength of the 
staggered magnetic order. The solid lines depicting phase 
boundaries may not be taken as true transitions at any fi- 
nite doping: impurity-induced AF order within a spin liq- 
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uid crosses smoothly to impurity-suppressed order within 
an AF, with no vanishing order parameter. Impurities 
may be considered as damaging both to the AF ordered 
phase, and the quantum spin hquid phase. 
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FIG. 5. Phase diagram for the 1/5-depleted, square-lattice 
antiferromagnet with doping by non-magnetic impurities. 

Experimentally, the situation in CaV409 may be more 
complex than that simulated here. Inelastic neutron scat- 
tering measurements , and subsequent refinements by 



comparison with simulations 1 13|, show that the magnetic 
system is composed of metaplaquettes with J2 — J = 
14meV, which interact with their neighbors via mutually 
frustrating couplings Ji = J[ = 0.49J, and with next 
neighbors by a further weak coupling — 0.25 J. Be- 
cause of the frustrated nature of the subdominant inter- 
actions Ji, it would appear reasonable that the effective 
coupling J' between metaplaquettes is small, and thus 
that the physical system is well in the PRVB regime; 
this deduction is consistent with the robust spin gap, 
A = 9.4meV, of CaV409. Doping with non-magnetic im- 
purities could be effected by random replacement of some 
spin-1/2 V^"*" ions by non-magnetic Ti^+. Our simula- 
tions of the uniform susceptibility and staggered magne- 
tization suggest that the weak AF order induced on dop- 
ing should be detectable in 5-10% Ti-doped CaV409, by 
measurements of x(2^)i Knight shift or magnetic Bragg 
scattering. 

In summary, we have analyzed the behavior on doping 
by non-magnetic impurities of the spin-liquid and ordered 
AF phases of the Heisenberg model on a 1/5-depleted, 
square lattice of spins S = 1/2. In the spin liquids, intro- 
duction of impurities creates effectively free spins which 
combine through exponentially weak interactions to cre- 
ate a long-ranged AF order at T = 0, coexisting with the 
gapped, spin-liquid phase. In the intrinsically AF regime, 
impurities cause a progressive destruction of the long- 
range order until it vanishes at the percolation threshold. 
The SSE QMC algorithm which we employ proves to be 
a powerful and efficient means of simulating the prop- 
erties of a spin system; these attributes are particularly 
important when many realizations of random impurity 
distributions are required. Finally, the unique magnetic 
plaquette structure of CaV409 provides a well-defined 
starting point for studying the effect of non- magnetic im- 



purities in spin liquids, and provides valuable insight for 
other 2d RVB compounds, such as frustrated Heisenberg 
systems, where a recombination of pruned bonds into new 
RVB clusters is also expected to occur. 
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